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We present a Generalized Langevin Equation for the dynamics of interacting semiflexible polymer 
chains, undergoing slow cooperative dynamics. The calculated Gaussian intermolecular center-of- 
mass and monomer potentials, wich enter the GLE, are in quantitative agreement with computer 
simulation data. The experimentally observed, short-time subdiffusive regime of the polymer mean- 
square displacements, emerges here from the competition between the intramolecular and the inter- 
molecular mean-force potentials. 

PACS numbers:61.25.Hq 



I 

o 
o 



> 
m 

O 



O 



X3 
O 

o 



X 



Despite its intrinsic complexity, the dynamics of low 
molecular weight unentangled polymer fluids has been 
considered for a long time to be a well understood prob- 
lem in polymer physics. It is commonly accepted that 
when the degree of polymerization, N, does not exceed 
the entanglement value, N^, a molecule is free to dif- 
fuse in the liquid, and follows Fickian dynamics. In this 
case, the Rouse approach successfully describes several 
key features of the polymer dynamics. These include 
the scaling with N of the bulk viscosity and the single- 
chain diffusion coefficient Duouse — kBT{(^N)~^ , with 
kB the Boltzmann constant, T the temperature, and 
the monomer friction coefficient . 

Upon closer examination, however, the dynamics of 
unentangled polymer fluids still presents important unre- 
solved questions. In the Rouse theory, local intermolecu- 
lar interactions are completely neglected; the single chain 
dynamics is driven by intramolecular entropic restoring 
forces and segmental friction, while the dynamics of the 
surrounding chains provide a heat bath. In this descrip- 
tion the polymer center of mass is free to diffuse following 
Brownian dynamics, and the mean-square displacement 
is expected to scale linearly in time at all time scales. 
In reality, a single polymer in a fluid spans a volume 
V oc Rg, with Rg = \/N 1/^/6 the polymer radius of gy- 
ration, and I the statistical segment length. Inside this 
volume are contained an average of n oc ^/N chains, that 
interact with each-other through the potential of mean 
force . The range of the potential is the same as that 
of the correlation hole Q in the pair correlation func- 
tion g{r), a distance of order Rg in polymer fluids [Q. 
Such interactions on length scales shorter than or equal 
to Rg leads to the Rouse equation's failure to adequately 
describe short-time polymer dynamics. 

The inconsistency of the Rouse equation with com- 
puter simulation data of short-time unentangled 
and entangled polymer dynamics in the melt state, 
and in concentrated solutions 0J^, has been known for 
many years. In the short-time regime, t <^ T^ouse = 
2R'i/{n^D Rouse), the single-chain center-of-mass (c.o.m.) 



mean-square-displacement, Ai?^(i), of unentangled poly- 
mers shows anomalous diffusion. For times t <C Tuouse 
AR^ (t) crosses over from the short-time ballistic dynam- 
ics to long-time Fickian diffusion (Ai?^(i) oc t) through 
an intermediate regime (Ai?2(t) cx iJ|,i- Only 

at long times {t > Tuouse) does the system obey Rouse 
dynamics and Fickian diffusion is recovered. Entangled 
polymer fluids show the same behavior in the short- 
time regime where Rouse dynamics is expected to hold 
The short-time anomalous dynamics only appears 
in melts and in solutions above the polymer overlap 
concentrations 0,0,^, further supporting the hypothe- 
sis that their short-time anomalous dynamics is a conse- 
quence of intermolecular interactions. 

In a recent paper, we derived a Generalized Langevin 
Equation for the Cooperative Dynamics (CDGLE) of in- 
teracting, completely flexible polymer fluids where 
explicitly intramolecular and intermolecular forces are 
included. The CDGLE is here implemented to treat 
the heterogeneous dynamics of unentangled semiflexible 
polymer melts with flnite-size chains. The theory is found 
to be in excellent agreement with computer simulation 
data performed by G. Crest and coworkers 

The CLE for a group of n molecules undergoing slow 
cooperative dynamics is derived from the flrst-principle 
Liouville equation by projecting the dynamics of the 
entire fluid onto the phase-space of the slow variables, 
i.e. the coordinates of n molecules statistically com- 
prised inside the volume defined by the range of the 
mean-force potential, n = pRgP = 0{^fN), with p the 
monomer density. An effective segment a in molecule i is 
driven by the intramolecular potential, — [3~^lri^ {r'^^\t)) 
with r«(i) = {v^{t),vf{t),..v^{t),..v^{i)}, the in- 
termolecular time-dependent potential of mean-force, 
— /3^^/ng[r(^-' (t), r^*^) (t)], and the random interactions 
with the surrounding fluid, where the projected random 
force is given by ¥^'^^\t). The segment position ra\t) 
follows the equation of motion : 
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n9(r^-''^W,r('='W)]+F?«(i) , 



(1) 



with 13 = {kBT)-\ 

The effective monomer friction coefficient, Ce//i is 
given by the linear combination of the bare friction, 
Co = /? < • Fa > /3, conventionally assumed to 
be the Rouse friction, and the intramolecular (i) and in- 
termolecular (j i) memory functions containing time 
and space correlation of the random forces: 
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+ ?EE / dt<F«(0).F?(^)(t)> . (2) 

The coupling of the random forces acting on different 
molecules is discarded in the traditional single-chain ap- 
proaches 0, and it is negligible for systems of uncorre- 
lated molecules, for which the pair distribution function 
g{r) « 1. 

With the purpose of describing the dynamics of finite- 
size semiflexible polymers, Eq.( 0) is reformulated in a 
matrix formalism, and semifiexibility is introduced fol- 
lowing the approach developed by Bixon and Zwanzig 
at the simple Rouse level [0. The new equation of 
motion applies to the overdamped dynamics of a fluid 
of interacting molecules of any kind, including semiflexi- 
ble polymers of finite length. It recovers the traditional 
single- molecule equation [0 in the limit n = 1, and/or 
with the fiuid described as a continuum, g(r,t) — 1. In 
the case n = 2 and = 1, it recovers the diffusive 
equation for binary interacting particles | |T^ . Thus our 
equation correctly reproduces the dynamics of systems 
characterized by dominant intramolecular or intcrmolcc- 
ular interactions. 

Since Eq.( ^ is obtained by projecting the fluid dy- 
namics onto the extended set of intra- and intcrmolcc- 
ular slow variables, the memory function contributions 
are minimized and can be neglected in first approxima- 
tion 1^. In addition, if only the intramolecular poly- 
mer coordinates are chosen as slow variables, and the 
memory function is discarded, this procedure recovers 
the traditional Rouse equation. The memory function 
contribution to the intramolecular equation becomes rel- 
evant for melts of high molecular weight polymers, where 
it describes the cross-over from unentangled to entangled 
polymer dynamics [^ . 

The GLE is decoupled in center-of-mass diffusion and 
intramolecular dynamics by a transformation of variables 
from beads to bond monomer coordinates Ip3 . The solu- 



tion of the c.o.m. diffusive equation requires the deriva- 
tion of an analytical expression for the c.o.m. intermolec- 
ular potential. We start from the four-point distribution 
function which relates a pair of monomers belonging to 
different polymers, and their c.o.m. positions. The dis- 
tribution is factorized in the product of two-point distri- 
bution functions and integrated over the monomer coor- 
dinates g{R) « / dva'^iri'^) J drb^{ri'''>)g{ra,b). In the 
special case of polymer molecules we adopt the thread- 
PRISM expression of the monomer pair distribution func- 
tion m , and the Gaussian intramolecular distribution , 
^'(ra) = [3/(27ri?2)]3/2ea;p[-3r2/(2i?2)]. The effective 
potential acting between the centers-of-mass of a pair of 
molecules inside a region of slow cooperative dynamics is 
approximated by 



w{R) = - In g(i?) « f ^ 1. 

[1 - 108^-2(^*)-2(^)-l]g-3flV(4i?.^,) 



(3) 



w{R) is Gaussian and finite for all distances, with a range 
on the order of the radius of gyration. At full polymer- 
polymer overlap, w{0) decreases with increasing reduced 
fiuid density, p* — pl^, reflecting the transition towards 
an incompressible system. w{Q) also decreases with in- 
creasing polymer molecular weight due to the increasing 
interpenetrability of the polymer chains. Figure 1 shows 
a comparison between Eq. ( ^ and the potential extracted 
from the trajectories of molecular dynamics computer 
simulations for melts of unentangled polyethylene chains 
The simulations are performed at constant volume 
and constant temperature, using the experimental den- 
sities at atmospheric pressure. We analyze trajectories 
from unentangled polyethylene samples with increasing 
degree of polymerization, as reported in Table I (the en- 
tanglement degree of polymerization Ne = 136). The 
only fltting parameter in the quantitative comparison, is 
the potential at complete polymer overlap, w'(0), which 
is compared in Table I to w(Q). All the samples are in 
good agreement with Eq.( ||). The value of w'(0) agrees 
quantitatively with the simulations for the high N sam- 
ples. In the low N samples the presence of the fine struc- 
ture at short distance, due to the local monomer packing, 
doesn't allow us to make a quantitative comparison. 

The data show a decrease of w(0) with increasing 
molecular weigth or density, in a greement with Eq.( |^), 
and with a recent study on dilute and semidilute poly- 
mer solutions [T^ . Eq.( |3|) represents an improvement 
on the mean-field formula derived by Flory and Krig- 
baum, which predicts the opposite trend [ p^ . However 
in solution both simulations and renormalization group 
calculations shows that the iV — *■ oo scaling regime is fi- 
nite. In our study of polymer melt dynamics, the small 
range of N available does not allow us to investigate this 
feature. 
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From Eq.( we derive the intermolecular force be- 
tween the c.o.m. of two polymers. In Eq.( ^ this 
is moltiphed by the number of interacting polymers 
G{t)R{t) = -VNp*dW{t)/dR{t). The time evolution 
of the intermolecular center-of-mass potential follows the 
time-dependence of the interpolymer distance. With the 
purpose of obtaining an analytical solution of Eq.( |^), we 
approximate B?(t) with its statistical average over the 
polymer configuration space, 

Gm(t) « - 1087r-2(p*)-2(7V)-i] 

e-^<R\t)>/(^K)R[t) . (4) 

When Eq.(l) is introduced in Eq. ( |]) , the latter reduces 
to a set of coupled Rouse equations with explicit inter- 
molecular contributions 

^^dR|^ = GW[RS„.(t) - W] + F\t) , (5) 

which obey the fluctuation-dissipation condition < 
Fl^it) ■ Fi^{t') 6Cof3-^Sa.jS,,6{t - t'), where a,7 = 

x,y,z. Rt(t) — n^^J2k=i'^'^m{t) is the coordinate of 
the c.o.m. of the dynamical aggregate. 

The advantage of adopting the approximate form of 
the intermolecular force is that the GLE becomes solv- 
able using the standard techniques of transformation to 
normal mode coordinates. The set of coupled equations, 
Eq.( reduces to two GLEs in the relative and coop- 
erative variables, from which the relative many-chain in- 
termolecular distance < R^{t) nAR^{t) — QDm.c.t, 
where AR^(t) is the single chain mean-square displace- 
ment. The many-chain cooperative diffusion coefficient, 
Dm.c. = {f3p*Ny^Coy^ = D Rouse /ip*VN), approaches 
zero in melts of high molecular weight polymer chains, 
at low temperature, and/or at high liquid density. Since 
the intermolecular force depends at each instant on the 
value of the intermolecular distance < R^{t) >, we 
solve Eq.( ^ self-consistently for small time increments 
{At = 10^^ ps). By this procedure we take into account 
the change in the intermolecular interaction with the dy- 
namical evolution of the system. 

Using Eq.( we calculate the c.o.m. mean-square dis- 
placement, AR'^{t) =< {Rc.m.{t) -i?c.m.(0))^ >. Since 
the appearance of anomalous c.o.m. dynamics is related 
to the amplitude of G{t) (when G{t) — ^ single chain 
uncorrelated dynamics is recovered) the anomalies are 
confined in the time scale t < Tuouse, which is the time 
necessary for the system to travel outside the range of 
the mean- force potential > R^)- 

We test the validity of our approach by comparing the 
c.o.m. mean-square displacement calculated from the so- 
lution of Eq.( 1^), with the data from computer simula- 
tions for the samples previously analyzed. The values of 
p, T, iV, / and Co are the input to the equation together 



with the effective value of the potential at initial time 
and zero intermolecular distance, ^'(O), calculated from 
the fitting of the potential. Here the statistical segment 
is much larger than the chemical bond because of the 
intrinsic stiffness of the polyethylene molecule. Figures 
2 and 3 show that the diffusive equation with the soft- 
core intermolecular potential reproduces quantitatively 
the anomalous subdiffusive c.o.m. dynamics for all of 
the samples investigated, and in the complete range of 
time scales. The slope of the subdiffusive c.o.m. regime 
is reported in Table I, and it is found to decrease with in- 
creasing degree of polymerization. The degree of dynam- 
ical anomaly increases with increasing N and decreasing 
T, as the system moves towards its glass transition. For 
the range of N investigated in this study a scaling regime 
has yet to be reached. 

To describe the polymer flexibility we adopt a freely ro- 

o 

fating chain model for segments of length I = 1.53 A with 
an effective stiffness parameter g = 0.74 {g —< 1; • l^+i > 
/P, li = rj+i — r.i). The value of g is obtained from 
the polymer mean-square end-to-end distance calculated 
from the simulation data, and reproduces the experimen- 
tal characteristic ratio for a finite polyethylene chain. 
The introduction of the explicit description of the local 
stiffness does not affect the dynamics at the c.o.m. level. 
This is because both the flexible model with renormalized 
segments and the the explicit semiflexible description are 
devised to reproduce the global chain properties. In fig- 
ures 2 and 3 the c.o.m. diffusion is identical if we use a 
semiflexible or a renormalized flexible model. 

The monomer dynamics depends strongly on the local 
chemical structure, and chain stiffness. However, at the 
monomer level, the intermolecular contributions are neg- 
ligible. Our mode analysis shows that only the lowest 
index, global modes are affected by the intermolecular 
potential in agreement with previous studies (|,@,|l^. In 
analogy with the static picture, we argue that the in- 
sensitivity of the local dynamics to the intermolecular 
forces is due to the presence of similar and compensat- 
ing monomer-monomer intramolecular and intermolecu- 
lar excluded- volume interactions In the regime of 
t < TRouse , the data show a dynamics slightly faster than 
the Rouse prediction of Ar^(t) oc t^^^ . This effect is pre- 
dicted by our diffusion equation with included local poly- 
mer semiflexibility. Figure 3 shows that the monomer 
mean-square displacement calculated with the semiflexi- 
ble approach is in better agreement with the simulations 
than the Rouse function. The inclusion, or exclusion, 
of the intermolecular contribution, does not modify the 
monomer diffusion. 

In conclusion we show that the center-of-mass inter- 
molecular potential in polymer melts can be described at 
a good level of approximation by an effective Gaussian 
potential. The introduction of the local intermolecular 
forces in the GLE strongly improves the description of 
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the global dynamical properties of polymer melts in the 
short-time domain. The intermolecular forces, however, 
do not affect the dynamics on the local monomer scale, 
which is instead strongly dependent on the local polymer 
semiflexibility. 
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TABLE CAPTION 

Table I: Simulation parameters, and fitting parame- 
ters. 

FIGURE CAPTIONS. 

FIG. 1: Comparison between the analytical expres- 
sion of the soft-core center-of-mass Gaussian potential, 
Eq.( H), and computer simulations data, vs. the c.o.m. 
interpolymer distance normalized by the polymer end-to- 
end distance, Rete — VnI- 

FIG. 2: Center-of-mass mean-square displacement as 
a function of time. Best fit of the molecular dynamics 
simulation data (filled circle) with the Rouse equation 
(dashed line), and with the intermolecular diffusion equa- 
tion, Eq.( 1^), for polymer melts (full line). The short- 
dashed lines indicate the longest Rouse relaxation time, 
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FIG. 3: Monomer and center-of-mass mean-square dis- 
placements as a function of time. Best fit of the molecu- 
lar dynamics simulation data (filled circle for c.o.m. and 
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open squares for monomer dynamics) with the Rouse 
equation (dashed hne), and with the diffusion equations 
for semiflexible polymer melts (full line). The short- 
dashed lines indicate the longest Rouse relaxation time, 
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